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I .  INTRODUCTION 


This  report  discusses  research  on  the  development  of  a  statistically 
based  model  for  predicting  human  head/neck  Impact  acceleration  Injury. 

Model  development  is  focused  on  those  situations  in  which  the  torso  is  well- 
restrained,  but  the  head  and  neck  are  unrestrained.  One  example,  of  course, 
is  the  Navy  pilot  in  the  cockpit  of  his  plane. 

The  objective  of  the  research  discussed  here  is  the  development  of  a 
model  that  can  adequately  predict  the  probability  of  head/neck  injury  based 
on  head  dynamic  response  data.  Once  this  objective  is  met,  the  model,  and  the 
information  it  provides,  should  be  a  major  component  in  development  of  im¬ 
proved  restraint  systems  and  other  protection  methods. 

In  order  to  provide  a  basis  for  comparison,  three  classes  of  prediction 
models  are  considered  in  this  report.  One  class  is  based  on  head  dynamic  re¬ 
sponse  variables  only,  another  is  based  on  sled  acceleration  profile  terms  only, 
and  the  third  class  is  based  on  the  combined  set  of  independent  variables.  All 
of  the  prediction  models  are  of  the  same  functional  type  as  those  considered 
in  previous  Desmatlcs  technical  reports  [6],  [7],  [8]: 

-1 

P(x)  -  {1  +  exp[-(gQ  +  EgjX^)]  } 

where : 

X  *  (x^, . . . ,Xj^)  denotes  the  set  of  independent  variables  considered, 

(8^,8^, . . .  ,Bj^)  denotes  a  set  of  parameter  values, 
and  P(x)  denotes  the  true  probability  of  injury  corresponding  to  x. 

This  type  of  model  was  previously  applied  to  observed  data  from  a  set  of 
28  -G^  accelerator  runs  involving  subhuman  primates  (Rhesus  monkeys)  with 
securely  restrained  torso  and  unrestrained  head  [7].  The  data  was  collected 
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by  the  Naval  Blodynamlcs  Laboratory  (NBDL)  as  part  of  its  research  effort  on 
Impact  acceleration  injury  prevention.  The  NBDL  data  base  now  consists  of  68 
-G^  accelerator  runs.  (The  28  runs  that  were  analyzed  previously  constitute 
a  subset  of  the  existing  data.)  In  addition  to  examining  this  larger  data 
set,  this  report  considers  additional  independent  variables.  The  variables 
comprising  the  former  data  set  and  the  additions  to  the  new  one  are  listed  in 
Table  1. 

The  data  base  is  used  to  develop  the  "best"  one-variable,  two-variable 
and  three-variable  models  for  each  of  the  three  classes.  In  the  context  here, 
the  "best"  model  Is  the  one  which  maximizes  the  log  likelihood  function  at 
each  stage  [4].  (In  all  cases,  the  contribution  of  additional  terms  beyond 
the  three-variable  models  was  negligible.)  The  predictions  from  the  "best" 
models  are  compared  with  the  observed  results  to  evaluate  performance. 


r-.. 


Sled  Profile  Variables 

*1.  Peak  acceleration  (G) 

*2.  Rate  of  acceleration  onset  (G/sec) 

3.  Duration  of  peak  acceleration  (msec) 


Head  Dynamic  Response  Variables 

2 

*1.  Peak  resultant  angular  acceleration  (rad/sec  ) 

2 

*2.  Peak  resultant  linear  acceleration  (m/sec  ) 

*3.  Peak  resultant  angular  velocity  (rad/sec) 

2 

4.  Peak  x-component  of  angular  acceleration  (rad/sec  ) 

5.  Peak  y-component  of  angular  acceleration  (rad/sec  ) 

6.  Peak  z-component  of  angular  acceleration  (rad/sec^) 

2 

7.  Peak  x-component  of  linear  acceleration  (m/sec  ) 

2 

8.  Peak  y-component  of  linear  acceleration  (m/sec  ) 

2 

9.  Peak  z-component  of  linear  acceleration  (m/sec  ) 

10.  Peak  x-component  of  angular  velocity  (rad/sec) 

11.  Peak  y-component  of  angular  velocity  (rad/sec) 

12.  Peak  z-component  of  angular  velocity  (rad/sec) 


*Denotes  a  variable  in  the  former  data  set. 


Table  1:  Independent  Variables  Available  for  Model  Building 
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II.  MODEL  CONSTRUCTION 


As  mentioned  in  the  previous  section,  the  current  NBDL  data  base  con¬ 
sists  of  68  observations.  However,  because  of  missing  data  on  five  of  these 
runs,  only  63  observations  were  used  in  model  construction.  Since  some  of  the 
monkeys  were  run  more  than  once,  dependence  exists  in  the  data.  However,  in 
model  development  the  assumption  is  made  that  any  resulting  bias  in  the  param¬ 
eter  estimates  is  small.  In  fact,  if  there  is  a  bias,  it  should  result  in  a 
model  that  overpredicts  probabilities.  This  is,  of  course,  the  best  direction 
for  model  bias,  since  it  provides  an  extra  margin  of  safety. 

Since  the  occurrence  or  nonoccurrence  of  injury  is  difficult  to  determine, 
the  criterion  of  fatality  is  used  in  the  model  building  process.  The  models 
are  thus  fatality  prediction  models.  The  data  for  all  63  observations  is  pre¬ 
sented  in  Table  2.  In  this  table,  the  observed  probability  of  fatality  for  a 
given  accelerator  run  is  denoted  by  1  for  a  fatal  run  and  0  for  a  nonfatal  run 

A.  DETERMINATION  OF  APPROPRIATE  MODELS 


A  forward  selection  method  [1,2]  was  used  to  determine  the  inclusion  of 
important  variables.  Importance  of  each  of  the  variables  was  determined  by 
likelihood-ratio  tests  [3]  that  are  used  in  conjuction  with  nested  models. 
This  involved  computing  the  following  quantities: 

Lj  * -2  log  likelihood  for  model  containing  (x^ , . . .  ,Xj^) . 

and  L2  ■ -2  log  likelihood  for  model  containing  (Xj^ . ^k-fl^‘ 

Under  the  null  hypothesis  that  the  additional  variable  provides  no  im¬ 

provement  in  the  model,  the  statistic  L^^  -  L2  has  an  approximate  Chi-square 
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Table  2:  The  Data  Set 
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The  Data  Set  (continued) 


distribution  with  one  degree  of  freedom.  The  hypothesis  may  be  tested  by  com¬ 
paring  the  value  of  -  L2  with  the  upper  percentage  points  of  the  Chi-square 
distribution. 

B.  HEAD  DYNAMIC  REPSONSE  VARIABLES 

As  previously  mentioned,  up  to  three  head  dynamic  response  variables  were 
considered  in  model  development.  As  a  consequence  of  the  forward  selection  pro¬ 
cedure  used,  (i.e.,  at  each  stage  the  variable  that  maximized  the  log  likelihood 
function  was  entered  next)  the  variables  In  the  "best"  one-variable  and  two-vari¬ 
able  models  constituted  a  subset  of  the  variables  chosen  for  the  "best"  three- 
variable  model,  A  statistical  analysis  of  the  data  indicated  that  the  best  one- 
variable,  two-variable,  and  three-variable  models  are  those  based  on,  respective¬ 
ly,  the  three  sets 
(1) 

(2)  Xj^,  X2 

(3)  Xj^,  X2,  x^ 

where 

x^  denotes  the  peak  z-component  of  head  linear  acceleration  measured  in 
meters/sec^, 

x_  denotes  the  peak  head  resultant  linear  acceleration  measured  in  meters/ 

^  7 

sec^ , 

and  x^  denotes  the  peak  y-component  of  head  angular  acceleration  measured  in 
radians/sec^ . 

Because  of  the  nesting  in  these  models,  the  relative  contribution  of  each 

of  the  variables  may  be  tested.  The  log  likelihoods,  the  Chi-square  values  and 

associated  p-values  presented  in  Table  3  contain  the  relevant  information.  In 

the  first  stage,  Xj^  was  tested  to  determine  whether  it  significantly  improved  a 

model  which  assumed  constant  probability  over  all  the  values  of  the  three  head 
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Variable  Set 


2  Loe  Likelihood 


Chi-Square 


-value 


Constant  Only  55.13 


25.58  29.55  0.00 

Xj^,X2  24.06  1.52  0,22 

Xi,X2,X3  20.64  3.42  0.06 


denotes  peak  z-component  of  head  linear  acceleration 
y^  denotes  peak  head  linear  resultant  acceleration 
y^  denotes  peak  y-component  of  head  angular  acceleration 


Table  3:  Head  Dynamic  Response  Variable  Sets  with 
-2  Log  Likelihood  and  Chi-Square  Values 


-8- 


dynamic  response  variables.  The  observed  Chi-square  value  of  29.55,  which  is 
statistically  significant  at  the  0.001  level,  indicated  that  this  variable  did 
result  in  an  improved  model. 

The  second  stage  of  testing  involved  consideration  of  the  addition  of 
another  variable  to  the  model  which  included  only  variable  x^.  Variable 
was  the  next  to  enter  the  model.  The  addition  of  X2,  which  resulted  in  an 
observed  Chi-square  value  of  1.52  and  an  associated  p-value  of  0.22,  did 
little  to  Improve  the  model.  However,  when  x^  was  added  to  the  model  con¬ 
taining  x^and  X2,  the  Chi-square  value  of  3.42  and  its  corresponding  p-value 
of  0.06  indicated  that  there  was  an  enhancement  to  the  model. 

Thus,  based  on  the  data  available,  the  best  one-variable,  two-variable  and 
three-variable  head  dynamic  response  models  are: 

?(x^)  =  {l  +  exp[-(-5.7852  -  .0048099x2)]}"^  (1) 

P(Xj^,X2)  =  {l  +  exp[-(-6.4795  -  .0035446x^+  .0008661x2)])“^  (2) 

P(Xj^,X2,x^)  =  {l  +  exp[-(-8.1485  -  .004019x^+  .001901x2  +  .000117x2)]}"^ 

(3) 

where  P(x)  denotes  the  predicted  probability.  The  discussion  of  the  contri¬ 
bution  supplied  by  each  of  the  variables  to  the  model  indicates  that  either 
model  (1)  or  (3)  could  be  chosen  for  prediction  purposes.  Table  4  presents, 
for  both  models,  a  comparison  of  observed  (i.e.,  0  or  1)  and  predicted  prob¬ 
ability,  where  the  observations  are  arranged  in  order  of  Increasing  predicted 
probability  for  model  (1) . 
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0.0012 

-165.0 

580.0 

2750.0 

LX1898 

A04101 

0 

0.0076 

0.0014 

-190.0 

670.0 

3800.0 

LX1083 

A03921 

0 

0.0078 
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1500.0 
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0 
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2700.0 
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0 
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0 
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-290.0 

625.0 
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0 
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0 
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0 
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0 
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8200.0 
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0 

0.0434 
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-560.0 

1630.0 
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LX 1900 

A04101 

0 

0.0499 

0.0177 

-590.0 

1540.0 

10000.0 

LX1901 

A04101 

c 

0.0748 

0.0327 

-680.0 

1680.0 

10000.0 

LX3702 

AR8872 

0 

0.0748 

0.0163 

-680.0 

1060.0 

6000.0 

LX1891 

A03943 

0 

0.0818 

0.0312 
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1950.0 

15500.0 

LX3693 

AR0012 

0 

0.0818 

0.0390 

-700.0 

1500.0 

6200.0 

LX3701 

AR8872 

0 

0.0893 

0.0237 

-720.0 

1130.0 

5250.0 

LX3715 

ARf:AC2 

0 

0.0893 

0.0297 

-720.0 

1340.0 

6700.0 

LX3714 

ARrAP2 

0 

0.1109 

0.0888 

-770.0 

1660.0 

3700.0 

LX3192 

AR8R66 

I 

0. 1133 

0.2301 

-775.0 

2380.0 

6000.0 

LX  3694 

AR0012 

0 

0.1259 

0.0659 

-800.0 

1500.0 

4900.0 

LX1359 
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0 

0.1369 

0.2149 
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2760.0 
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LX1363 
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1 
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0.6059 

-840.0 

2000.0 

-12000.0 

LX3188 

AR8863 

1 

0.1549 

0.1606 

-850.0 

2600.0 

16000.0 

LX 13 62 

A03935 

0 

0.1613 

0.3750 

-860.0 

4100.0 

31000.0 

LX1902 

A041O1 

0 

0,1747 

0.3401 

-880.0 

1770.0 

-5000.0 

LX3706 

AR8802 

0 

0. 1965 

0.0402 

-910.0 

1460.0 

12500.0 

LX1903 
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0 

0.2203 

0.  5167 

-940.0 

1 oon.O 

-5600.0 

LX3705 

ARR802 

0 

0.2287 

0.0345 

-050.0 

1390.0 

1 4100.0 

LX37no 

AP.8790 

\ 

0.3137 

0.2043 

-1040. 0 

2200.0 

13500.0 
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AR8857 

0 

0.3347 

0.4493 
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0 
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-1350.0 

2400.0 

21000.0 
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A03933 

0 

0.6700 

0.5932 

-1  350.0 

2550.0 

15000.0 

LX1905 
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1 

0.8068 

0.9851 

-1 500.0 

2810.0 

-8300.0 

LX1365 

AO  3921 

1 

0.8416 

1 . 0000 

-1550.0 

9200.0 

22600.0 

LX1895 

A03951 

1 

0.8958 

0.8845 

-1650.0 

2880.0 

16500.0 

LX1896 

AO 3 94 6 

1 

0.9162 

0.6763 

-1 700.0 

2490.0 

23000.0 

LX1897 

AC3896 

1 

0.9162 

0.9062 

-1700.0 

3050.0 

inooo.o 

LX1360 

AC4099 

I 

0.9974 

n.  9009 

-2440.0 

52  00  .0 

1.9500.0 

Table  A:  A  Comparison  of  Observed  and  Predicted 
Probabilities  for  Head  Dynamic  Response 
Models  (1)  Xj  and  (3)  ,  x^. 
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C.  SLED  PROFILE  VARIABLES 


The  same  technique  used  in  the  previous  section  was  also  employed  for 
choosing  the  "best"  one-variable,  two-variable,  and  three-variable  sled 
profile  models.  The  sled  profile  variables  under  consideration  here  are 
denoted  by  z^t  and  z^,  where: 

z^  is  the  peak  sled  acceleration  measured  in  G's, 

Z2  is  the  duration  of  peak  measured  in  milliseconds, 
and  z^  is  the  rate  of  sled  acceleration  onset  measured  in  G/sec. 

As  Table  5  Indicates,  the  "best"  one-variable  model  is  based  on  z^^.  The 
Chi-square  value  of  38.43  for  and  its  corresponding  p-value  of  0.001  re¬ 
veals  that  peak  sled  acceleration  is  essential  to  the  sled  profile  model. 

The  resulting  model  is  given  by: 

P(z^)  -  {1  + exp (-(-12.10+  .11462z^)]}~^ 

The  second  stage  of  testing  provided  the  following  "best"  two-variable 
model : 

P(Zj^,Z2)  -  {l+cxp[-(-6.3859+  .11239Zj^  -  .30402Z2)  ]  }~^  (5) 

However,  the  addition  of  Z2,  duration  of  peak,  did  little  for  the  betterment 

of  the  model,  with  a  Chi-square  value  of  1.37  and  a  p-value  of  0.24. 

Subsequently,  in  the  third  stage,  z^  was  added  to  the  model  containing  z^ 

and  Z2-  The  resulting  three-variable  model  is: 

P(z^,z^,Zj)  =  {1 +  exp(-(-5.6780+  .10391Zj^  -  .34906z2  +  .00003557z2)  ] 

(6) 

The  Chi-square  value  of  0.01  for  z^  is  evidence  that  this  variable  is  not  an 
Important  addition  to  the  model. 

As  indicated  by  these  results,  the  "best"  sled  profile  model  is  the  one- 

variable  model  (4)  containing  only  peak  sled  acceleration.  Table  6  presents, 
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Variable  Set  -2  Log  Likelihood  Chi-Square 

Constant  Only  55.13 

16.70  38.43 

15.33  1.37 

Zj.z^.Zj  15.32  0.01 


denotes  peak  sled  acceleration 
Z2  denotes  duration  of  peak 
Zj  denotes  rate  of  onset 


Table  5:  Sled  Acceleration  Profile  Variable  Sets  with 
-2  Log  Likelihood  and  Chi-Square  Values 


-value 


0.00 


0.24 


0.93 
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0 

0.0754 

83.7 

LX1891 

A03943 

0 

0.0762 

83 . 8 

LX3709 

AR879n 

\ 

0.1166 

87  .  n 

UX31 88 

493863 

\ 

0.4695 

1  04 . 5 

LX31 92 

AR8866 

\ 

0.4924 

105.3 

LX1362 

A03935 

0 

0. 4''R1 

105.5 

LX1359 

A04099 

0 

0,5381 

106 . 9 

LX1894 

An3'’33 

0 

0.5861 

108.6 

LXl 365 

A03921 

1 

0 . 5888 

108.7 

LX1R93 

A03924 

0 

0,6351 

110.4 

LXl 363 

A03935 

1 

0,8806 

123.0 

LX1905 

4041  01 

\ 

0.9159 

126.4 

LXl 360 

404099 

\ 

0.9305 

128.2 

LX1895 

A03951 

] 

0.9469 

130.7 

LXl 896 

A03946 

\ 

0.9508 

131.4 

LXl 897 

A03896 

\ 

1.0000 

192.9 

Table  6;  A  Comparison  of  Observed  and  Predicted  Probability 
for  Model  Based  on  Peak  Sled  Acceleration 


-13- 


for  this  model,  a  comparison  of  observed  (l.e.,  0  or  1)  and  predicted  prob¬ 
ability,  where  the  observations  are  arranged  in  order  of  increasing  predicted 
probability. 


D.  COMBINED  HEAD  AND  SLED  VARIABLES 


The  complete  set  of  head  dynamic  response  and  sled  profile  variables 
was  also  used  in  the  development  of  a  prediction  model.  As  Table  7  indicates, 
y^  =  (peak  sled  acceleration)  provided  the  best  fitting  one-variable  model, 
which  is  given  by  (A) .  The  next  term  to  enter  the  model  was  y^  (the  peak 
z-component  of  head  angular  velocity) .  This  variable  had  a  Chi-square  value 
of  2.42  with  an  associated  p-value  of  0.12.  The  resulting  model  is: 

P(yj^,y2)  -  {l+exp[-(-13.582+  .12863yj-  .0082432y2)  ] (7) 

The  "best"  three-variable  model  was  obtained  by  introducing  y^  H  Z2 
(duration  of  peak  sled  acceleration)  into  the  model  containing  y^  and  y2. 

The  resulting  model  is: 

P(yj^.y2.y3)  -  {l  +  exp[-(-6.9958+  .12234y^-  .008002y2  -  .37337y3)  ]  }"^  (8) 
However,  the  Chi-square  value  of  1.51  and  its  corresponding  p-value  of  0.22 
indicate  that  duration  of  peak  did  not  improve  much  on  the  "best"  two-variable 
model. 

Thus,  the  "best"  combined  (head  dynamic  response  and  sled  profile)  model 
appears  to  be  the  two-variable  model  containing  peak  sled  acceleration  and  the 
peak  z-component  of  head  angular  velocity.  Table  8  shows  the  agreement  between 
predictions  and  observations  that  is  obtained  for  this  model. 
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Variable  Set 


-2  Log  Likelihood 


Chi-Square 


-value 


Constant  Only  55.13 

16.70  38.43  0.00 

y^,y2  14.28  2.42  0.12 

12.77  1.51  0.22 


y^  denotes  peak  sled  acceleration 

y^  denotes  ppak  z-component  of  head  angular  velocity 
y^  denotes  duration  of  peak  sled  acceleration 


Table  7:  Combined  Variable  Sets  with  -2  Log  Likelihood 
and  Chi-Square  Values 
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74 . 7 

-28.0 

LXl  <>03 

A04101 

u 

0.0232 

75.3 

-19.0 

LX19nn 

A041 01 

0 

0.0323 

/4 . 8 

-68.0 

LX3186 

LXSIO" 

AR8857 

ARR790 

0 

1 

0.2150 

0.2194 

82  .  7 

87 . 9 

-200.0 

-122.0 

LX1342 

A03935 

0 

0.2854 

105.5 

110.0 

LX3188 

ARRS63 

1 

0.2929 

104.5 

90.0 

LX1804 

A03933 

u 

0.3540 

108.6 

120.0 

LXl 359 

A04099 

0 

0.4194 

106.9 

60 . 0 

LX3192 

AR8S66 

1 

0.6414 

105.3 

-75.0 

LXl 893 

A03024 

n 

0.8214 

110.4 

-110.0 

LX1895 

A0395I 

1 

0. 8295 

130.7 

2  00 . 0 

LX1363 

A03935 

1 

0.8747 

123.0 

36.0 

LyiS96 

AO  3 04 6 

1 

0.9432 

131.4 

62 . 0 

LXl 905 

A04ini 

\ 

0.9522 

126.4 

-38.0 

LX1365 

A03921 

A04O99 

1 

0.9639 

108.7 

-350.0 

LXl 360 

I 

0.9921 

I  28.2 

-233.0 

LXl 897 

A03896 

1 

1 . 0000 

192.9 

20.0 

Table  8:  A  Comparison  of  Observed  and  Predicted  Probability 
for  Model  Based  on  Peak  Sled  Acceleration  and  Peak 
z-Component  of  Head  Angular  Velocity 
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E.  CLASSIFICATION  OF  OBSERVATIONS 


The  predicted  probabilities  from  the  "best"  models  developed  In  the 
previous  sections  can  be  used  to  classify  observations  Into  groups.  In  other 
words,  an  observation  can  be  classified  as  nonfatal  If  the  predicted  probability 
Is  less  than  or  equal  to  some  specified  cut-off  value.  In  particular,  each  such 
value  yields  a  classification  matrix  of  the  form  given  In  Figure  1.  From  this 
figure,  the  following  probabilities  can  be  defined: 

■  Prob (observe  fatality  |  predict  nonfatallty)  =B/(B  +  D) 

and 

Prob(observe  fatality  |  predict  fatality)  = A/ (A + C) . 

Ideally,  It  Is  desired  to  have  P^^  =  0  and  P2“l.  Of  course,  Pj^  Is  the  more 
critical  probability  of  the  two. 

Graphs  of  P^^  versus  P2  (as  a  function  of  the  cut-off  value)  compare  the 
performance  of  the  "best"  one-varlable ,  two-varlable,  and  three-variable  models. 
For  example.  Figure  2  compares  the  models  based  on  the  head  dynamic  response 
variables  with  those  based  on  sled  acceleration  variables.  The  Improvement 
that  Is  obtained  by  going  from  the  "best"  one-varlable  to  the  "best"  three- 
variable  model  (In  terms  of  approaching  the  Ideal  situation,  l.e.,  P^ > 0  and  P2 
■  1)  can  be  seen  graphically  within  each  variable  set. 
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V  Predicted 

Observed 

Fatality 

Nonfatality 

Fatality 

A 

B 

Nonfatality 

1 

C 

D 

-  B/<B+D) 
P^-A/CA+C) 


Figure  1:  Classification  Matrix 


HI  denotes  best  one-variable  head  model 
H3  denotes  best  three-var table  head  model 
SI  denotes  best  one-variable  sled  model 
S3  denotes  best  three-variable  sled  model 

Figure  2:  Comparison  of  Models 
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III.  SUMMARY 


Using  an  Identical  data  base,  three  different  models  were  constructed, 
one  based  on  sled  profile  variables,  another  based  on  head  dynamic  response 
variables,  and  the  last  one  comprised  of  the  combined  set  of  Independent 
variables.  The  "best"  head  dynamic  response  model  appeared  to  be  the  three- 
variable  model  containing  the  peak  z-component  of  head  linear  acceleration, 
peak  head  linear  resultant  acceleration,  and  the  peak  y-component  of  head 
angular  acceleration.  The  "best"  sled  acceleration  profile  model  was  the  one- 
varlable  model  consisting  of  peak  sled  acceleration  alone.  The  "best"  combined 
(head  dynamic  and  sled  profile)  model  was  the  two-varlable  model  consisting  of 
peak  sled  acceleration  and  the  peak  z-component  of  head  angular  velocity. 

The  statistical  technique  of  testing  the  contribution  of  successive  terms 
In  nested  models  [3]  cannot  be  employed  for  models  Involving  different  variables 
(l.e.,  no  formal  test  exists  for  determining  whether  or  not  one  model  provides 
a  significant  Improvement  over  another).  However,  a  relative  assessment  of  the 
various  models  can  be  made  on  the  basis  of  the  log  likelihood  values.  In  par¬ 
ticular,  for  models  containing  the  same  number  of  variables,  the  one  which  max¬ 
imizes  the  log  likelihood  value  would  be  favored.  In  this  regard.  It  can  be 
seen  from  Table  9  that  the  three-variable  head  dynamic  response  model  does  not 
do  any  better  than  the  one-varlable  model  based  on  peak  sled  acceleration  alone. 
Similarly,  It  appears  that  the  combined  two-varlable  model  does  better  than  the 
three-variable  head  dynamic  response  model.  In  addition,  graphs  of  the  prob¬ 
ability  of  correct  classification  (l.e.,  observing  fatality  given  that  fatal¬ 
ity  Is  predicted)  Indicate  that  this  probability  Is  maximized  sooner  for  the 
combined  two-varlable  model. 

There  still  remains  the  question  as  to  why  the  head  dynamic  response  models 
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S.  Variable 
^\Type 

Best 

Models 

Head 

Dynamic 

Response 

Sled 

Acceleration 

Profile 

Combined 
Var iable 
Set 

One- Variable 

-12.79 

-8.35 

-8.35 

Two-Variable 

-12.03 

-7.66 

-7.14 

Three- Var lab le 

-10.32 

-7.66 

-6.38 

Table  9:  Log  Likelihood  Values  for  Best  Models 
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did  not  perform  as  well  as  the  models  involving  sled  profile  variables.  It 
was  speculated  in  [7]  that  this  result  may  be  due  to  any  or  all  of  the  follow¬ 
ing:  (a)  the  wrong  variables  were  being  extracted  from  the  head  dynamic  re¬ 
sponse  time  traces,  (b)  inaccurate  measurements  were  being  made  on  the  correct 
variables,  and  (c)  the  small  sample  size  had  produced  a  spurious  result. 

Since  the  sample  size  appears  to  be  sufficiently  large  here,  it  is  be¬ 
lieved  that  (c)  can  be  ruled  out.  In  addition,  (b)  also  appears  to  be  an  un¬ 
likely  explanation,  since  the  effect  of  minor  measurement  inaccuracies  would 
most  likely  be  negligible  in  larger  samples.  However,  it  is  still  possible 
that  there  exists  more  valuable  Information  that  can  be  extracted  (via  the 
method  of  principal  components  [5],  for  example)  from  the  head  dynamic  re¬ 
sponse  time  traces. 

This  latter  contention  is  borne  out  by  the  fact  that  the  head  dynamic  vari¬ 
able  found  to  be  most  important  was  the  peak  z-component  of  head  linear  ac¬ 
celeration,  which  was  not  available  for  consideration  as  a  variable  in  the 
original  report  on  model  development  [7].  This  indicates  that,  in  some  sense, 
the  current  set  of  twelve  head  dynamic  response  variables  provides  better  in¬ 
formation  than  the  original  set  of  three  variables. 
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